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Abstract Cladistic analysis of the morphological relationship of 40 populations of Paeonia sect. 
Moutan DC. (including all wild species) was carried out with an attempt to have a better under- 
standing of the phylogeny of tree peonies. Twenty-five morphological characters were selected and 
used in unweighted pair-group method using arithmetic averages (UPGMA), the neighbor-joining 
(NJ) and the maximum parsimony (MP) tree construction executed with PAUP" program. The 
topologies of the trees obtained were largely in congruence and topological differences occurred only 
between MP tree and the other two in a clade comprising five species, P. suffruticosa, P. jisha- 
nensis, P. qiui, P. rockii and P. ostii, which were considered as closely related species based 
on both morphological and cytological characters. In both NJ and UPGMA trees, all the populations 
of each species formed a monophyletic clade, and the clades of P. ludlowii, P. qiui and P. de- 
composita. had strong bootstrap supports ( > 9096) , the clades of P. rockii, P. delavayi had mod- 
erate bootstrap supports ( > 7096) , and the clades of P. suffruticosa, P. ostii and P. jishanensis 
had weak bootstrap supports (51% — 7096) . In MP tree, however, the corresponding bootstrap val- 
ues of all the clades mentioned above were generally smaller than those in NJ and UPGMA trees and 
the populations of P. delavayi, P. ostii and P. suffruticosa failed to form monophyletic clades. 
Based on the above results, interspecific relationships of wild tree peonies were discussed and a phy- 
logeny of sect. Moutan was reconstructed for the first time. 

Key words Paeonia, Paeonia section Moutan , morphology, population, phylogeny. 


Tree peonies belong to sect. Moutan DC. in the genus Paeonia L. and are all endemic to 
China (Pan, 1995). In the past decade, important progresses have been made on the taxonomy of 
tree peonies (Hong, 1997a, 1997b, 1998; Hong D-Y et al. , 1996, 1998a, 1998b; Hong & Pan, 
1999a, 1999b; Hong T et al., 1992; Pei & Hong, 1995). However, little work was done on the 
phylogeny of sect. Moutan (Zou et al., 1999). Pei and Hong (1995) considered that P. qiui Y. 
L. Pei & D. Y. Hong is morphologically most closely related to P. jishanensis T. Hong & W. Z. 
Zhao. Hong and his co-workers (Hong et al., 1998a) stated that P. delavayi Franch. is most 
closely related to P. ludlowii D. Y. Hong, and the cultivated tree peony (P. suffruticosa Andr.) 
may have been derived from P. suffruticosa ssp. yinpingmudan D. Y. Hong, K. Y. Pan & Z. 
W. Xie (Hong et al., 1998b). Using a combined data matrix of both external morphological and 
pollen characters, Yuan and Wang (1999) built a UPGMA tree of five wild tree peony species, 
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i.e. P. decomposita Hand.-Mazz., P. ostii T. Hong & J. X. Zhang, P. rockii T. Hong & J. 
J. Li, P. qiui and P. jishanensis, and a hybrid, P. x yananensis T. Hong & M. R. Li. Be- 
cause no bootstrap values were given in their tree, the relationships among the species obviously 
need to be corroborated. Using RAPD markers of 59 populations from seven wild species, P. de- 
composita, P. ostii, P. rockii, P. qiui, P. jishanensis, P. delavayi and P. ludlowii, and 
one hybrid, P. x yananensis , Zou and her co-workers (Zou et al., 1999) built a UPGMA tree of 
all the species studied, but the tree produced has a topology not in accordance with that of Yuan and 
Wang (1999). Recently some very important progresses have been made in studies on phylogeny 
and hybrid speciation in sect. Paeonia using DNA sequences (Sang, 2002; Sang et al., 1995, 
1997a, 1997b; Ferguson & Sang, 2001; Tank & Sang, 2001). However, information on interspe- 
cific relationships within sect. Moutan provided by those studies is still quite limited because of the 
limited sampling (Sang et al., 1995, 1997b; Ferguson & Sang, 2001; Tank & Sang, 2001). 
Therefore, the main purpose of the present study is to analyze morphological relationship among the 
species of sect. Moutan using different phylogenetic analysis methods based on population sam- 


pling, with an attempt to have a better understanding of the phylogeny of wild tree peonies. 
1 Materials and Methods 


1.1 Materials 

The specimens used in this study were all collected from natural populations, except for P. 
suffruticosa ssp. suffruticosa which was a photo of the type stored in the Herbarium, Institute of 
Botany, the Chinese Academy of Sciences, Beijing (PE). Forty populations (each with three or 
more specimens) representing all eight species and their corresponding subspecies (Hong & Pan, 
19993) were investigated. In addition, two specimens respectively from the two herbaceous sections 
of the genus Paeonia, P. brownii Douglas ex Hooker and P. lactiflora Pall., were used as func- 
tional outgroups in tree constructions. The origin of the materials used in the study is presented in 
Table 1. 

1.2 Methods 

The names of species used in this study follow the classification of Paeonia sect. Moutan made 
by Hong and Pan (1999a). Twenty-five morphological characters were selected and used in this 
analysis for their importance in species delimitation. The characters, their states and code numbers 
are given in detail in Table 2. NEXUS Data Editor ( Version 0.5.0) was used to record raw data 
and create data matrix for tree-building calculations. 

The maximum parsimony (MP), neighbor-joining (NJ), and unweighted pair-group method 
using arithmetic average (UPGMA) trees were constructed using PAUP” 4.0b4a computer program 
(Swofford, 2000). During tree construction process, all characters were unordered and equally 
weighted, with multi-states interpreted as uncertainty and “?” treated as missing information. Parsi- 


mony analysis was performed by heuristic search with tree-bisection-reconnection (TBR) of branch 
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Table 1 Information of populations of species in Paeonia sect. Moutan used in this morphological analysis 

No. Names of populations " Origin Voucher ? Determined by 
1 lacti Weichang, Hebei Chengde Team 28 D. Y. Hong 
2 brownii USA - D. Y. Hong 
3 suffruticosa Type photo - - 

4 yinpingmudan- HNSX Song Xian, Henan H97010 D. Y. Hong 
5 yinpingmudan-AHCH Chaohu, Anhui 9701 D. Y. Hong 
6 jishanensis-SX ]S Jishan, Shanxi 96053, PB89501 D. Y. Hong 
7 Jishanensis- HNJY Jiyuan, Henan H97001, 940220 D. Y. Hong 
8 Jishanensis-ShXYA Yan' an, Shaanxi 94008 D. Y. Hong 
9 decomposita-KD Kangding, Sichuan 95080A D. Y. Hong 
10 decomposita-MEK Barkam, Sichuan PB85025 D. Y. Hong 
ll decomposita-DB Danba, Sichuan 95050-3B D. Y. Hong 
12 decomposita- JC Jinchuan, Sichuan 9113003 D. Y. Hong 
13 rotundiloba-MX Mao Xian, Sichuan PB89701, 95015-2B D. Y. Hong 
14 rotundiloba-LX Li Xian, Sichuan 9112004 D. Y. Hong 
15 rotundiloba- WC Wenchuan, Sichuan 9111009 D. Y. Hong 
16 rotundiloba- HS Heshui, Sichuan 95017-4-9 D. Y. Hong 
17 ostii-ShXHL Huangling, Shaanxi H97069 D. Y. Hong 
18 ostii-HBBK Baokang, Hubei 97052A D. Y. Hong 
19 ostii-ShXTBS Taibaishan, Shaanxi PB85052 D. Y. Hong 
20 ostii-HNNX Neixiang, Henan H97021A D. Y. Hong 
21 qiui-HBBK Baokang, Hubei H97027 D. Y. Hong 
22 qiui-HBSX Songbai, Hubei PB85034 D. Y. Hong 
23 qiui- HNXX Xixia, Henan PB88305 D. Y. Hong 
24 taibaishanica-ShXMX Mei Xian, Shaanxi H97058 D. Y. Hong 
25 taibaishanica-GSTS Tianshui, Gansu 13 D. Y. Hong 
26 taibaishanica-XZ ( C) Xizang (cultivated) H96051 D. Y. Hong 
21 rockü-HBBK Baokang, Hubei H97051 D. Y. Hong 
28 rockii-HBSBX Songbai, Hubei PB88029 D. Y. Hong 
29 rockii-HNNX Neixiang, Henan H97016-1B D. Y. Hong 
30 rockii- HNSX Song Xian, Henan 94002 D. Y. Hong 
31 ludlowii-MLI Mainling, Xizang H96182 D. Y. Hong 
32 ludlowii-ML2 Mainling, Xizang H96014 D. Y. Hong 
33 ludlowii-ML3 Mainling, Xizang H96030 D. Y. Hong 
34 delavayi-YNZD Zhongdian, Yunnan H97128 D. Y. Hong 
35 delavayi- YNL] Lijang, Yunnan H97095-11 D. Y. Hong 
36 delavayi-YNKM Kunming, Yunnan H97077-4 D. Y. Hong 
37 delavayi-SCDW Dawu, Sichuan 95063 D. Y. Hong 
38 delavayi- SCY] Yajiang, Sichuan 95070-3 D. Y. Hong 
39 delavayi-SCDK Panzhihua ( Dukou) , Sichuan 11458 D. Y. Hong 
40 delavayi-SCYY Yanyuan, Sichuan 97110-3 D. Y. Hong 
4l delavayi-XZBM Bomi, Xizang 379 D. Y. Hong 
42 delavayi-XZLZ Nyingchi, Xizang 3229 D. Y. Hong 


1) The names of populations are epithets of the species (subspecies) followed by the abbreviation of the population origin. For exam- 
ple, "yinpingmudan- HNSX" means the population of P. suffruticosa ssp. yinpingmudan. from Song Xian, Henan Province. 2) All 


the vouchers are preserved in PE. 
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Table 2 Morphological characters, their states and code numbers used in data matrix 








No. Characters Character states and their code numbers 
1 leaf type tri- or tetra-compound (0); bi- or tri-pinnate (1); biternate (2) 
2 leaflet number more than 35 (0); more than 10 (1); 5-9 (2) 
3 shape of leaflets rounded to elliptic (0); ovate to long-ovate (1); ovate-lanceolate to 
lanceolate (2) 
4 terminal leaflets entire or rarely lobed (0); cleft less than 1/2 of leaflet length (1); cleft 
more than 1/2 of leaflet length (2); divided to the base of leaflets (3) 
5 lateral leaflets entire or rarely lobed (0); cleft less than 1/2 of leaflet length (1); cleft 
more than 1/2 of leaflet length (2); divided to the base of leaflets (3) 
6 leaflet apex rounded to obtuse (0); acute (1); acuminate (2) 
7 leaflet base decurrent (0); cuneate (1); broadly-cuneate to rounded (2) ; 
8 leaf indumentum beneath glabrous (0); hairy only along veins (1); obviously hairy throughout the 
lower surface (2) 
9 flower number > 1 (0); solitary (1) 
10 flowering status more or less pendent (0) ; erect (1) 
11 bract number 1 (0); 2-3 (1); 4-5 (2) 
12 petal number single (0); double (1) 
13 petal color white (0); yellow (1); pink to red (2); purple (3) 
14 spots at petal base absent (0); red or pink-purple (1); dark purple (2) 
15 filament color white (0); yellow (1); pink to red purple (2) 
16 carpel number >5 (0); 2-4 (1); 1 or rarely 2 (2) 
17 carpel indumentum glabrous (0); denseły hairy (1) 
18 stigma color pale yellow to yellow (0); red (1); purple (2) 
19 flower disc leathery (0); fleshy (1) 
20 degree of envelopment of carpels by dise whole (0); only half of carpels (1); only the base of carpels (2) 
21 disc color white (0); pale yellow to yellow (1); red to purple (2) 
22 follicle shape ellipsoid (0); cylindrical (1); long ovoid (2) 
23 root shape ordinary (0); tuberous (1) 
24 reproduction obligately sexual (0); facultatively vegetative (1) 
25 habit subshrubs (0); herbs (1) 


swapping algorithm, the MULTREES, and Accelerated Transformation ( ACCTRAN) options. 
UPGMA and NJ trees were constructed using distance measure of mean character difference. The 
number of 1000 replicates is used for all bootstrap tests. 


2 Results 


2.1 Topology of trees and bootstrap support of species clades 

The NJ, UPGMA and one of the MP trees are shown in Figs. 1 — 3. The topologies of the three 
trees were generally in congruence, and the incongruences only occurred between MP tree and the 
remaining two in the clade consisting of five closely related species, i.e. P. suffruticosa, P. 
jishanensis, P. qiui, P. rockii and P. ostii (Figs. 1 - 3). In both the NJ and the UPGMA 
trees, all the populations of each species formed a monophyletic group with bootstrap supports at 
above 50% . On these two trees, the clades of P. ludlowii, P. giui and P . decomposita were 
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Table 3 The data matrix of 42 populations and 25 morphological characters used in tree construction analysis 





No. Names of populations! Code numbers of character states 

1 lacti 22(012) (23) (23)2010120(03) (01)1(01)001210101 
2 brownii 22(01) (23) (23) (02)0001103110021210101 
3 suffruticosa 221 (12) (01) 110112130201 100207 00 

4 yinpingmudan-HNSX 221(12)121011203020110020? 00 

5 yinpingmudan-AHCH 221 (12) 12101 1200020120020? 00 

6 jishanensis-SXJS 221 (123) 1(12) (12)211000201 100201110 

7 jishanensis-HNJY 22(01)111(12)211100021 110020110 

8 jishanensis-ShXYA 22(01)211(12)2111032(12)0110020110 

9 decomposita-KD 00(01) (12) (01)21011203000010100000 

10 decomposita-MEK 000(23) (12)21011(12)0(23)0(01)0000100000 
11 decomposita- DB 00(01)3(12)21011203000010100000 

12 decomposita- JC 0003111011(12)03000010100000 

13 rotundiloba-MX 00(01) (13) (01)1101110300101020(01)000 
14 rotundiloba-LX 00(01) (12) 111011103001010100000 

15 rotundiloba- WC 00(01)1111011103001010100000 

16 rotundiloba-HS 000(23) 1 (12) 101110300101010(01 )000 

17 ostii- HBBK 112002101 10000201 10020000 

18 ostii-ShXHL 2(1 + 2)2002101120(0 + 2)020110020000 
19 ostii-ShXTBS 22(02) (12)021011200021 110020000 

20 ostii- HNNX 222102101 12000201 10020000 

21 qiui-HBBK 221101201110(03) 120120022? 10 

22 qiui-HBSX 221(01)012111103120120022? 10 

23 qiui-HNXX 2210012111103120120022? 10 

24 taibaishanica-ShXMX 110(12) (12)21211200200100010100 

25 taibaishanica-GSTS 1101021211200200100010100 

26 taibaishanica-XZ ( C) 11011(12)1111200201100010100 

27 rockiü-HBBK 112202111110(03)210100010000 

28 rockii- HBSBX 11200211 1 11002001001 10000 

29 rockii- HNNX 11200211 1 1100200100010000 

30 rockii- HNSX 112002121 1100200100010000 

31 ludlowii-MLI 2212220000101012001211100 

32 ludlowii-ML2 2212220000101012001211100 

33 ludlowii-ML3 2212220000101012001211100 

34 delavayi-YNZD 221(23) (23)200(01)010(012)1(12)10212(12)0110 
35 delavayi- YNLJ 221222001010(23)121011220110 

36 delavayi-YNKM 221332000010(13)121021220110 

37 delavayi-SCDW 2223320000(12)0(12) (01)210112(12)2110 
38 delavayi-SCYJ 2223320010201 (01 )2102122(02)110 

39 delavayi-SCDK 222332001020( 123 ) (01)210212(12)2110 
40 delavayi-SCYY 2223320000102221021220110 

4l delavayi-XZBM 222(23) (23)20000(12)0102101122(02)110 
42 delavayi-XZLZ 222(23) (23)20000( 12)01021011227 110 


1) The names of populations ars defined the same as in Table 1. 
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Fig. 1. The unweighted pair-group method using arithmetic averages (UPGMA) tree of the tree peony populations studied in the 
present analyses. Numbers above branches indicate the mean distance; numbers below branches are bootstrap values. 


strongly supported (bootstrap values > 90% ) , and those of P. rockii and P. delavayi were moder- 


ately supported (bootstrap values > 70% ) (Figs. 1, 2). However, in the MP tree, the bootstrap 
figures of all the species clades were generally smaller than those in the NJ and UPGMA trees and 
the populations of P. delavayi, P. ostii and P. suffruticosa failed to form their own monophyletic 


groups (Fig. 3). 


2.2 Phylogenetic relationship among species of Paeonia sect. Moutan DC. 
In the trees produced (Figs. 1 — 3), P. delavayi and P. ludlowii first branched out from the 
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Fig. 2. The neighbor-joining tree of the tree peony populations studied in the present analyses. 
Numbers above branches indicate the mean distance; numbers below branches are bootstrap values. 


other species within sect. Moutan. Because the clade was strongly supported by bootstrap values 
(96% in the NJ and UPGMA trees, 90% in the MP tree), a closer relationship between P. lud- 
lowii and P. delavayi than between either of them and any other species in sect. Moutan was 
strongly supported. The populations of the two subspecies of P. decomposita first formed their own 
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Fig. 3. The maximum parsimony tree of the tree peony populations studied in the present analysis. Numbers above branches 
indicate the number of character changes that occurred along each branch; numbers below branches are bootstrap values. 


monophyletic group respectively (NJ and UPGMA trees), and then together formed a monophyletic 
clade with a strong bootstrap support (9796 in both NJ and UPGMA trees, 9096 in MP tree), 
which in turn consisted of a sister relationship with the clade comprising five species, viz. , P. suf- 
fruticosa, P. jishanensis , P. qiui, P. rockii and P. ostii. For this species group, their inter- 
specific relationships were not well resolved by this analysis because of the topological differences 
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among the trees obtained and the lack of robust bootstrap support for the branches (Figs. 1-3). 
The populations of the two subspecies of P. rockii formed their own monophyletic groups ( bootstrap 
> 7896), and they then formed a monophyletic clade (bootstrap > 70% ) which was sister to the 
other four species in both the NJ and UPGMA trees (Figs. 1, 2). This sister relationship, howev- 
er, was not topologically supported by MP tree (Fig. 3). In both the NJ and UPGMA trees, the 
populations of P. jishanensis and P. qiui first formed two independent clades, and then joined in 
a monophyletic group with 60% bootstrap value (Figs. 1, 2), but this group collapsed again in the 
MP tree (Fig. 3). In addition, the populations of P. rockii and P. ostii formed a monophyletic 
group, which was then joined by P. suffruticosa in MP tree (Fig. 3). In the same way, these rela- 
tionships were not supported by both the NJ and UPGMA trees. The relationship among these tree 
peony species obviously needs further investigation. 

3 Discussion 

Following de Candolle (1824), Stem (1946) divided the genus Paeonia into three sections, 
sect. Paeonia, sect. Onaepia, and sect. Moutan. The sect. Moutan , with four species and two 
varieties, was further differentiated into two subsections, viz., subsect. Vaginatae and subsect. 
Delavayanae (Stern, 1946). Based on a series of field investigations and taxonomic studies on wild 
tree peonies (Hong, 1997a, 1997b, 1998; Hong et al., 1998a, 1998b; Hong & Pan, 1999b) , 
Hong and Pan (1999a) suggested a new classification system for sect. Moutan, in which eight 
species (with three of them each having two subspecies) were recognized. The present study showed 
that all the populations of each species in Hong and Pan’s (1999a) system formed their own mono- 
phyletic clades in all three trees (Figs. 1-3), except for the populations of P. delavayi and P. 
suffruticosa in the MP tree, and they together formed a monophyletic group having a sister relation- 
ship with the outgroups. Therefore, this result is in consistent with the taxonomic treatment of 
species in sect. Moutan made by Hong and Pan (19992) . In addition, the subdivision of P. rockii 
into two subspecies by Hong (1997a) is also supported in this analysis. 

Up to date, our understanding on the phylogeny of sect. Moutan was still limited and sporadic 
(Sang et al., 1995, 1997b; Zou et al., 1999; Yuan & Wang, 1999; Ferguson & Sang, 2001; 
Tank & Sang, 2001). A close relationship between P. ludlowii and P. delavayi and their rela- 
tionship with all the other species in sect. Moutan are well recognized by both morphological study 
(Hong & Pan, 19992) and molecular evidence (Zou et al. , 1999) , and are also confirmed by the 
present analysis. Hong and Pan (1999a) considered that three groups could be distinguished in 
sect. Moutan, with the first group comprising P. ludlowii and P. delavayi, the second P. de- 
composita , and the third P. suffruticosa, P. jishanensis, P. qiui, P. rockii and P. ostii. This 
division implies that P. decomposita has an intermediate position among the species of sect. 
Moutan . From Figs. 1 — 3 we can see that the suggested systematic position of P. decomposita is 
well supported and this result is also consistent with earlier conclusions (Tank & Sang, 2001; Yuan 
& Wang, 1999) . 
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Previous studies on sect. Moutan showed that the interspecific relationships among P. suffru- 
ticosa, P. jishanensis, P. qiui, P. rockii and P. ostii are problematic ( Hong & Pan, 1999b; 
Yuan & Wang, 1999; Zou et al., 1999). These five species were considered to be closely related 
in both morphological and cytological aspects and were treated as a species complex ( P. suffruticosa 
complex) by Hong and Pan (1999b). All the three trees presented in this work (Figs. 1 - 3) 
showed that the five species formed a monophyletic clade with P. decomposita as its sister group. 
Pei and Hong (1995) suggested that P. qiui is most closely related to P. jishanensis and both of 
them have a close relationship with P. suffruti- P jishanensis 
cosa. A close relationship between P. giui and A 
P. jishanensis is supported by the UPGMA analy- Bgint 
sis of morphological and palynological data (Yuan | 
& Wang, 1999), but not by molecular evidence "Po 
(Zou et al. , 1999) and only weakly supported by 
our NJ tree (Fig. 2). In addition, P. rockii and 


P. ostii showed a closer relationship than those 


P ostii 


P rockii 


between either of them and any other species with- 
in P. suffruticosa complex in both NJ and MP 
trees ( Figs. 2, 3), whereas P. suffruticosa 
showed a closer relationship with P. ostii than 
between P. rockii and P. ostii in UPGMA tree 
(Fig. 1). Based on previous works and the pre- 


P. decomposita 


P. delavayi 





P. ludlowii 


: : Fig. 4. A suggested interspecific relationship in Paeonia 
sent study mentioned above, a phylogenetic rela- |." Moutan DC. The solid lined clades are well s spanned ui 


tionship among species in Paeonia sect. Moutan the present study and the dotted lined clades need further evi- 


i i in Fi nfirm. 
is suggested, as is shown in Fig. 4. dence to confirm 
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为 药 属 牡 丹 组 基于 形态 学 证 据 的 


系统 发 育 关 系 分 析 
"3 周志 钦 EFE REL 


!( 中 国 科学 院 植物 研究 所 系统 与 进化 植物 学 重点 实验 室 “北京 100093) 
《西南 农业 大 学 园艺 系 ”重庆 400716) 


摘要 ”对 芍药 属 牡丹 组 Paeonia L. sect. Moutan DC. (全 部 野生 种 )40 个 居 群 进行 了 基于 形态 学 证 据 的 系 
统 学 分 析 , 试 图 建立 组 内 种 间 的 系统 发 育 关 系 。 利 用 PAUP* 计算 机 程序 分 别 构建 了 建立 在 25 个 形态 学 
性 状 基 础 上 的 所 有 研究 类 群 的 距离 树 (UPGMA、NJ) 和 最 大 简约 树 (MP)。 所 得 树 的 拓扑 结构 基本 一 致 ， 
差异 只 发 生 在 距离 树 和 简约 树 之 间 , 在 由 形态 和 细胞 学 关系 都 很 近 的 5 个 种 (牡丹 P，suffruicosa、 矮 特 
Ft P. jishanensis 、 卵 叶 牡 丹 P，giwi、 紫 斑 牡 丹 P. rocki MAF P，ostii ) 构 成 的 分 支 内 部 。 在 距离 树 
(UPGMA 和 NJ) 中 ,每 个 种 的 所 有 居 群 都 首先 形成 1 个 单 系 类 群 ,然后 才 与 其 他 种 的 居 群 相聚 。 大 花 黄 牡 
F P. ludowi 、 卵 叶 牡 丹 和 四 川 牡丹 P. decomposita 3 个 种 的 单 系 类 群 得 到 大 于 90% 的 自 展 值 支持 ,紫斑 
牡丹 和 演 牡 丹 P. delavayi 2 个 种 的 单 系 类 群 得 到 大 于 70% 的 自 展 值 支持 ,牡丹 \ 凤 丹 和 矮 牡丹 3 个 种 的 
单 系 类 群 得 到 51% - 70% 的 自 展 值 的 支持 。 在 简约 树 (MP) 中 ,上 述 所 有 种 的 自 展 支 持 值 都 小 于 距离 树 
中 的 相应 数值 。 同 时 牡丹 、 汗 牡丹 和 凤 丹 的 全 部 群体 未 能 形成 一 个 单 系 群 。 根 据 本 文 的 研究 结果 , 分 
析 讨 论 并 首次 提出 了 野生 牡丹 的 系统 发 育 关系 。 
KR VAR; 牡丹 ; 形态 学 ; BHR; 系统 发 育 


